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We investigate T-matrix for bound and continuous-spectrum states in the discrete oscillator 
representation. The investigation is carried out for a model problem - the particle in the field 
of a central potential. A system of linear equations is derived to determine the coefficients 
of the T-matrix expansion in the oscillator functions. We selected four potentials (Gaussian, 
exponential, Yukawa, and square-well ones) to demonstrate peculiarities of the T-matrix and 
its dependence on the potential shape. We also study how the T-matrix expansion coefficients 
depend on the parameters of the oscillator basis such as the oscillator length and the number 
of basis functions involved in calculations. 
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1. Introduction 

We are going to consider the convergence of a wave 
function and the T-matrix expansion in the oscillator 
functions. This consideration will be restricted to a 
model problem of the particle in the field of a central 
potential. The analysis will be done within the ma¬ 
trix form of quantum mechanics, which involves an 
infinite set of oscillator functions to realize a discrete 
representation. This matrix form is well known as the 
algebraic version of the resonating group method or 
J-matrix method. The methods were formulated in 
mm and Ena, respectively. Now, they are widely 
used to describe nuclear, atomic, and molecular sys¬ 
tems. In Ref. jSj, one can find progress in resolving 
the internal problems of the method and the numer¬ 
ous applications to solving real physical problems in 
different branches of quantum physics. 

Usually, the expansion of a wave function in the 
oscillator or any other square-integrable basis is con¬ 
sidered within the J-matrix. The discrete form of the 
T-matrix has not been investigated yet. For instance, 
in Ref. [S], the oscillator basis was used to construct 
wave functions and to calculate the phase shift for a 


Gaussian potential. The T-matrix was also obtained 
in that paper, but only in the momentum space. So, 
we are going to fill in this gap. We consider the states 
of both discrete and continuous spectra. However, the 
main attention will be paid to the scattering states. 

We are going to demonstrate the convergence of 
the T-matrix and how the convergence depends on 
the shape of a potential, on the oscillator length of a 
basis and the energy of the state. For this aim, we se¬ 
lected four potentials, which mimic different physical 
systems. Note that, by solving model problems, one 
can reveal the interesting features and peculiarities 
of systems under consideration, which are observed in 
more complicated and realistic systems. For instance, 
it was discovered in Refs. [7115]. while studying simple 
model problems, that the J-matrix in some cases suf¬ 
fers of a slow convergence. This means that one needs 
to involve a very large basis of oscillator functions to 
achieve the desired precision of calculations. An effec¬ 
tive method was formulated in [HHin]. It allows one 
to reduce the set of oscillator functions by three-five 
times in order to obtain the phase shifts with higher 
precision. 

The paper is organized in the following way. In 
Section |^ we make all necessary definitions and de- 
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duce the equations for the T-matrix expansion coeffi¬ 
cients. We briefly consider main equations, which are 
used to describe quantum systems. These equations 
will be transformed from a continuous (coordinate 
or momentum) representation to a discrete, oscilla¬ 
tor representation. The analysis of the T-matrix in 
the discrete representation is presented in Section 
for four potentials (Gaussian, Yukawa, exponential, 
and square-well ones). 

2. Model Formulation 

We make use of the system of units by selecting the 
constant h?/m = 1. This leads to a renormalization 
of the potential energy operator 

F(r)^gF(r). 

In this representation, the kinetic energy operator is 
T = — in coordinate space and T = in the 

momentum space, and the energy is if = In this 

paper, we will consider central potentials. Thus, the 
orbital momentum L is a good quantum number. 

2.1. Basic equations 

To determine the spectrum of bound states and their 
wave functions or to determine the scattering param¬ 
eters and the corresponding functions for continuous 
spectrum states, one should solve the Schrodinger 
equation 

[H - (r) = 0 (1) 

or the Lippmann-Schwinger equation 
(x) = IpkL (x) - 

— j G{x,x)V x^dxx dx. (2) 

The latter can be written in the coordinate space 
{x = r) or in the momentum space {x = p). In the 
coordinate space, we have 


( 3 ) 

( 4 ) 

( 5 ) 


’’PkL (r) = \l -kjL (kr), 
G{r,r) =jL {kr<)nL {kr-^), 
V {x, x'j = S (r) 


(r< = min (r, r), r-y — max (r, r)) and in the momen¬ 
tum space 


V’fcL (p) =Hp-k), 
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( 6 ) 


-1 


G {P,P) = - ^k'^ +'ie] Hp-p), (J) 


y {P,P) = -PP J JL (pr) V (r) Jl {pr) r^dr. (8) 

0 

Note that the transition from the coordinate space 
to the momentum one is determined by the Fourier- 
Bessel integral 


'^kL (p) = \j‘^P j jL {pr) 4'fcL {x) r'^dr. 


( 9 ) 


There is another equation, which is also used to 
determine the spectrum and the wave functions of 
bound and scattering states. This is the Lippmann- 
Schwinger equation for the half-off shell transition T- 
matrix (see, e.g., fTTlfT^ l 


OO 

tL {p, k) = Vl {p, k) + J dpp 


2 {P,P)tL {p, k) 
E — 4^ -I- ie 


■ ( 10 ) 


There are several equivalent definitions of the T-mat¬ 
rix, in particular, those, which involve integrals with 
the potential energy operator and the wave function 


th {p, k) = ^ ^p J drr (pr) V (r) (r) (11) 

0 

in the coordinate space or, in the momentum space. 


tL (p, ^) = J dp^y {P,P) 4'fcL {p). 
0 


( 12 ) 


We present the following important relation, which 
connects the T-matrix and the wave function in the 
momentum space: 

tL {p, k) 


^kL {p) =6{p-k)- J 


\p^ — \kP 


(13) 


By calculating the T-matrix, we can easily construct 
the wave function of a system in the momentum 
space. 

It is well known that, in order to determine the 
spectrum and wave functions by solving the Schro¬ 
dinger equation, one needs to impose the adequate 
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boundary conditions, while the necessary boundary 
conditions are included in the Lippmann-Schwinger 
equation. 

In the present paper, we use the standing-wave rep¬ 
resentation, which means that the asymptotic part of 
wave functions is 


(^) = JL (kr) - tanSL^J^k ul (kr), 

where Sl is the phase shift, and (kr) and ul (kr) 
are Bessel and Neumann functions, respectively. Ho¬ 
wever, the T-matrix is usually determined with the 
wave function in the running-wave representation. In 
this representation. 


W = Ui ^ (kr) - (kr) 


Here, Sl = exjp{2i5i,'\ is the scattering matrix, and 


{kr) = ^\l^k [-riL {kr) ± IJl {kr)] 

are incoming ^ {kr)) and outgoing {kf)) 

waves. 

Note that the factor \ -k in the definition of the 

Y 'TT 

asymptotic part of a wave function for the continu¬ 
ous spectrum state is chosen to normalize the wave 
function 5'fcL by the condition 

{'i>kL\'i’kL) =5ik-k). (14) 

It is easy to show that the T-matrix con¬ 

structed in the running-wave representation is con¬ 
nected with the T-matrix in the standing-wave 

one by the simple relation 

y(W) ^ _-gi5z. ginJ^T^^^l 

Thus, we prefer to work with the real function . 


2.2. Discrete representation 

To transform the Schrodinger equation and the 
Lippmann-Schwinger equation for a wave function or 
T-matrix to the discrete representation, we use a full 
set of oscillator functions {4>„i ( 2 ^, 6 )} in the coordi¬ 
nate space and {^nL {Py b)} in the momentum space. 
The explicit form of the functions is as follows: 

{r,b) = 

(15) 

P = r/b, 
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(p, 6 ) = (p^), 

p = pb. 

Here, 

J^nL = 


(16) 


2r (n -t 1) 


r (n -k T -t- 3/2) ■ 

These functions obey the completeness relations 


^ {r, b) (r, b) =6{r-r), 

n—0 

00 

<^nL {P, b) (p, b)=5{p-p). 


(17) 

(18) 


n—0 


To transform any equation for a wave function or 
T-matrix to the oscillator representation, we will use 
the orthogonality of the basis functions and the com¬ 
pleteness relations. We can also use the fact that all 
quantities, which appear in Eqs. Q, ([^, and (10 1 , 
can be represented as 

00 

F{x) = '^^nL{x,b)Fr,L{b), (19a) 

n—0 

00 

O (x) = ^ {x, b) (n d rnj ^mL {x, b), (19b) 

n—0 

00 

6 {x, x) = ^ (x, b) (n d (x, b), (19c) 


n—Q 


where x stands for r or p. In Eqs. (191, it is tacitly as¬ 
sumed that the function F (x), local operator O (x), 
and nonlocal one O (x, x) obey all necessary condi¬ 
tions to be expandable in the oscillator basis. 

We start the transformation with the Schrodinger 
equation. It is easy to verify that the Schrodinger 
equation in the oscillator representation is 


00 

EK" 


m—0 


H 


-ES„ 


Cm=0, 


( 20 ) 


where {Cm} are the wave function expansion coeffi¬ 
cients 

00 

4'fcL {r) = ^ Cm^nL {r, b), 

m—0 

Cm = i^nLl'^kL)- 


( 21 ) 
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Now, we turn our attention to the Lippmann-Schwin- 
ger equation for the wave function. The integral equa¬ 
tion ([^ is transformed to a system of linear equations 


OO OO 

E 

m—O m—0 


Cn = - E E 1^1 ^ ^ 4 


Here, is the wave function of free motion 

(Bessel function) ^ in the oscillator representa¬ 
tion. Due to the peculiarities of oscillator functions 
and the oscillator Hamiltonian, the expansion coeffi¬ 
cients coincide with oscillator functions in the 
momentum space (see more details in 

= {^uLli’kL) = ^uL{p = k,b). (23) 


Thus, we have to deal with an infinite vector. 

Second, we can represent the T-matrix as a matrix 


tLip,k)= ^ ^nL{P,b)tnmib)^mLik,b), (27) 

n,m—0 

where the matrix elements tnm (b) are determined as 

OO 

tnm {b) = j dpp^^nL {p, b) (p, k) X 


X ^mL {k, b) k^dk. 


(28) 


To solve system (22), one needs to calculate the ma¬ 


trix elements of the potential energy operator 

OO 

V rrij = J drr‘^^nL{r,b)V {r)^rnLir,b) = 

= J dpp^^nL {p, b) V (p, k) (k, b) k'^dk 
0 

and the matrix elements of Green’s function for the 
free motion Hamiltonian between oscillator functions 

(n |G| m) = 

OO 

= J dpp^<^nL {p, h) 

0 

OO OO 

= J J drr‘^^nLir,b)G{r,r)^mL{r,b)r^dr. (24) 
0 0 

In Ref. [TJ , one can find the explicit form of the ma¬ 
trix elements (n|G|m) and the recurrence relations 
they satisfy. 

There are two different ways to present the T-mat¬ 
rix in the discrete (oscillator) form. First, we can use 
the expansion 


We note that the expansion coefficients for the T- 
matrix in both representations depend on the oscil¬ 
lator length. In the next section, we will study how 
strongly the T-matrix expansion coefficients depend 
on the oscillator length b. 


By projecting Eq. (101 on the oscillator basis, we 


obtain the sets of linear inhomogeneous equations for 
the vector = t„L(6, k), 


tn = Vnib,k) + ^ {n 

m,rh—0 


V 


or 


E 


m—O _ 


dnm ^ ^ 


m—O 


V 


m) (m jGj m) tm (29) 


tm — (b, k ), 

(30) 


. |G| to) 


E-^p"^ + ie 

f t f 1L 

^mL (P, b)= ^ / 


2 

tnm — ^nm “1“ / ^ 

V 


and for the matrix (b)' 

to) (to |G| to) tfi 


(31) 


m,m—0 


or 


“'nm / \ ^ 

rh—0 L m—O 


E 

m—C 

Here, 


-E 


V 


. |G| to) 


tmm — ^nm- 


(32) 


th (p, ^) = E b) tnL {b, k). 

n—0 

(25) 

It is obvious that the expansion 
{tnL (b, fc)} are determined as 

OO 

coefficients 

tnL {b,k)= y dpp^$„L (p, 6)tL (p,/c). 

(26) 


VnL (b, k) = -k J drG^nL (r, b) V (r) jl (kr) = 
0 

OO ^ 

/ OO 

dpp^^nL (P, b) Vl {p, fc) = ^ (n |r| to) . 

n m=0 


Let us consider the T-matrix expansion coefficients 


(26). By using the definition of the T-matrix (11), 


this equation can be rewritten as 

tnL (b.k) = (^nL V 'hfei) = 
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oo 

= J drr‘^^nL {r, h) V^kL 
0 


or 

oo 

tnL {b, fc) = 

m—0 



m) C„ 


(33) 


(34) 


Here, {Cm} are expansion coefficients of the wave 
function '^kL (r) (see Eq. (HI)). 


2.3. Square-well potential 

To check our numerical results, we need to consider 
a potential, which admits a simple expression for the 
T-matrix. There is only a restricted number of two- 
body problems, which can be solved analytically, and 
for which the T-matrix can be obtained in a closed 
analytic form. For instance, the T-matrix in Ref. m 
was constructed for a delta-shell potential. 

One of these cases is the square-well potential 


V{r) 


(Vo r <a, 
(0 r > a. 


This case has been considered in detail many times. 
See, for instance. Ref. [16], where the off-shell T-mat¬ 
rix was obtained. We also consider this potential in 
detail, assuming that Vq is negative, and the potential 
is attractive. The obtained results will be intensively 
used in Section [S] 

To determine the T-matrix for a square-well po¬ 
tential, we have to obtain the wave function. The 
wave function for the potential in the internal region 
(0 < r < a) reads 

4'feL (r) = (r) = AkL]j^kjL (kor), (35) 

where Ak is the normalization constant, and 

ko = y/2 (E + Vo), k = V^. (36) 


In the asymptotic region (r > a) for the continuous- 
spectrum states, we have 




k jh (/cr)-tan(5i 



k (kr). 


(37) 


By matching the internal (351 and asymptotic (371 
parts of the wave function and their first derivatives 


as well, we can determine the phase shift Sl and the 
constant AkL- 

The T-matrix for the square-well potential is fully 
determined by the internal part (35) of the wave func¬ 
tion 

tL=tL {p, k) = ^J^VoAkL jdrr'^jL (pr) Jl (kor). (38) 


It is easy to verify that, for the zero value of orbital 
momentum L = 0, 


2 k 

to {p,k) = - — VoAk X 
TT ko 

p cos (pa) sin (koa) — ko sin (pa) cos (fepa) 

and 


(39) 


Ak = k/ [fcsin (koa) sin (ka) + ko cos (koa) cos (ka)]. 


On the energy shell, we have 


t (fc, k) 


— — tan Jo- 

TT ko 


(40) 


Note that expression ( [M| represents the T-matrix 
not only for scattering states (E > 0), but also for 
bound state(s) (—Vg < E < 0). One has to calculate 
the energy of a bound state (or the momentum ko) 
and the corresponding normalization factor Ak- 


3. Analysis of Results 

The numerical analysis of the T-matrix will be car¬ 
ried out for the s-state (L = 0) only, where the inter¬ 
action is more stronger than in other orbital states 
(L > 0) and is not diminished by the centrifugal 
barrier. Four potentials are used to study properties 
of the T-matrix. They are Gaussian (G), exponential 
(E), Yukawa (Y), and square-well (SW) potentials: 


Vir) 

Vir) 

V(r) 

V(r) 


(d1' 

Voexp|-^|, (E) 


(Vo r <a, 
\o r > a, 


(SW) 


(41) 
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Fig. 1. Convergence of the ground-state energy for the Gauss 
and exponential potentials 

In this section, we use the nuclear units for energy 
(MeV) and length (fm), so that the constant /m = 
= 41.47 MeV-fm2. 

For all potentials, we take the radius of the poten¬ 
tial 0 = 1 fm, and the depth Vq = —85 MeV. With 
such a choice of parameters, we obtain one bound 
state for all potentials, whose energies are listed in 
Table. It was established by the numerical solution 
of the set of equations ( [20| with a maximal num¬ 
ber of oscillator functions (see the text bellow). One 
can see that we obtained a deeply bound state for the 
Yukawa potential, a weakly bound state for the Gaus¬ 
sian potential, and a moderately bound state for the 
exponential and square-well potentials. The results 
demonstrated in Table are obtained with N = 300 
functions and the oscillator length 5 = a = 1.0 fm. 

In what follows, we will use four values of the os¬ 
cillator radius b — 0.5, 1.0, 1.5, and 2.0 fm to study 


Energy of the bound state 


Potential 

G 

E 

Y 

SW 

E, MeV 

-3.564 

-7.006 

-26.744 

-9.388 


the dependence of the T-matrix expansion coefficients 
on the oscillator radius. Calculations are organized in 
the following way. First, we construct the N x N ma¬ 
trix of the Hamiltonian, where 1< V < 500. (We 
do not dwell here on the calculations of matrix ele¬ 
ments of the potential energy operator between os¬ 
cillator functions. We make use of the technique of 
generating functions, details of which can be found 
in Refs. HZIIIE]). Second, we calculate the eigenval¬ 
ues (spectrum) and the eigenfunctions of the matrix 
by using the discrete form of the Schrodinger equa¬ 
tion ( |20[ ). Then we obtain the T-matrix of bound and 
pseudobound states in the oscillator representation, 
by using Eq. (34l. Third, we solve the system of lin¬ 
ear equation (221, which determine the wave function 
and the phase shift of the scattering state with fixed 
energy E, and then we construct the T-matrix by us¬ 
ing Eq. (34l. Finally, on the fourth step, we solve the 
set of linear equations (29) and obtain directly the 
T-matrix for a fixed energy E. We make use of the 
third way to check the correctness of the calculations 
of the T-matrix in the fourth way. 


In this section for the sake of convenience, we will 
denote the T-matrix as t (p, E) explicitly indicating 
the energy of the discrete or continuous-spectrum 
state. 


3.1. Convergence 

First of all, we consider whether the oscillator basis is 
large enough to provide the convergent results for the 
bound state energy and the phase shift of continuous- 
spectrum states. In Fig. we show the dependence 
of the bound state energy E on the number of os¬ 
cillator functions N involved in calculations. These 
results are obtained for two potentials (Gaussian and 
exponential) and for four different values of oscilla¬ 
tor length b. One can see that one hundred functions 
give a very stable value of ground state energy for 
all oscillator lengths. The dependence of the phase 
shift on the number of oscillator functions N used 
in calculations is demonstrated in Fig. The phase 
shift is determined for the energy E = 5.0 MeV. The 
exact values of phase shift are calculated within the 
variable phase method muni. To obtain the stable 
phase shift independent of N, we need to use more 
oscillator functions than for bound state calculations. 
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Fig. 2. Phase shift as a function of the number of oscillator 
functions involved in calculations. Results are obtained for the 
Gauss and exponential potentials with the energy E = 5.0 



n 

Fig. 3. T-matrix for bound states calculated with the Yukawa 
and exponential potentials 


3.2. Ground state 


In this section, we consider the T-matrix of bound 
states. In Fig. we show the T-matrix for bound 
states with the Yukawa and exponential potentials. 
To demonstrate the rate of decreasing of the T-matrix 
expansion coefficients, we display the renormalized 
expansion coefficients = tn {b,E) /to {b,E). For 
the Yukawa potential, the bound state is a deeply 
bound one with its energy to be —26.74 MeV. This 
explains why the T-matrix for this potential decreases 
very rapidly, as n increases. The common feature of 
the T-matrix for the Yukawa and exponential poten¬ 
tials is that the larger the oscillator length b, the 
slower is the decreasing of the T-matrix expansion 
coefficients. The same tendency is observed for the 
Gaussian and square-well potentials. 


With the help of relation (25), we can easily con¬ 


struct the T-matrix of a bound state in the momen¬ 
tum space. In Fig. we display the T-matrix t (p, E) 
as a function of the momentum p for the bound state 
for four potentials. These results are obtained with 
the oscillator length 6 = 1 fm and with N = 300 



Fig. 4. T -matrix of the bound state as a function of the 

momentum p calculated with four potentials 

oscillator functions. However, by using other values 
of oscillator length, we obtain the same results. It is 
worth to note that the T-matrix for a deeply bound 
state (Yukawa potential) is very dispersed in the mo¬ 
mentum space. This reflects the fact that the behav¬ 
ior of the T-matrix in the momentum space is de- 
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0 50 100 150 200 


n 

Fig. 5. T-matrix expansion coefficients for continuous-spect¬ 
rum states obtained with the Gaussian potential 


termined by the wave function in the internal and 
asymptotic regions (see Eq. In some cases, the 

internal part of wave functions gives a larger contri¬ 
bution than the asymptotic part. The same is true for 
the square-well potential, where the T-matrix is fully 
determined by the internal part of the wave function 
(as was pointed out above (see Eq. (38l)). Contrary 
to the case of the Yukawa and square-well potentials, 
the T-matrix for the Gaussian and exponential poten¬ 
tials is represented by the small values of momentum 
p: 0 < p < 7 fm“^. 


3.3. Continuous-spectrum states 

We selected five values of energy {E = 1, 5, 10, 15, 
and 20 MeV) to study the dependence of the T-matrix 
expansion coefficients on the energy in the continu¬ 
ous spectrum. The energy range 0 < E < 20 MeV se- 



n 

Fig. 6. The T-matrix expansion coefficients as a function of 
n calculated with the exponential potential 

lected in our calculations represents the region, where 
the effect of the potential energy is stronger than in 
the high-energy region. 

As we are interested in the study of the rate of de¬ 
creasing of , we will display the renormalized expan¬ 
sion coefficients tn = tn/to. In Fig. we show the T- 
matrix expansion coefficients for the Gauss potential 
for four values of energy E = 1, 5, 10, and 15 MeV. 
One can see that tn drops to zero very rapidly, as 
n is increased. The smaller the oscillator length, the 
faster is the decreasing of the T-matrix expansion co¬ 
efficients for the Gauss potential. Note that the rate 
of decreasing of almost independent of the energy 
used in our calculations. 

By comparing Fig. |^with Figs. ilZl and[^ we can 
see that the fast decrease of expansion coefficients is 
observed only for the Gaussian potential, while the 
other potentials exhibit a much more slower decrease. 

Gontrary to the Gaussian and square-well poten¬ 
tials, the T-matrix expansion coefficients for the ex¬ 
ponential and Yukawa potentials are more strongly 
dependent on the energy. This is explicitly demon¬ 
strated in Fig. where the expansion coefficients for 
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Fig. 7. Expansion coeficients tnito, obtained with the 
Yukawa potential 

the T-matrix are displayed for a fixed oscillator length 
b and for five different values of energy E. 

The T-matrix as a function of the momentum p for 
the square-well potential is shown in Fig. These 
results are obtained with N = 300 oscillator func¬ 
tions. In Fig. we compare the calculated T-matrix 
with the exact one (marked as Ext) represented by 
Eq. (39l. One can see that the calculated T-matrix 


P < 


y/AN -t 2L -k 3 


(42) 


Thus, the smaller the oscillator length, the larger 
is the range of momenta, which can be covered by 


Fig. 8. Square-well potential. The renormalized T-matrix 
expansion coefficients as a function of n 


N oscillator functions. Indeed, with b = 2.0 fm, the 
nonzero values of the T-matrix are obtained for p < 
< 17 fm“^, which is in accordance with formula (421. 


coincides with the exact one in a wide range of mo¬ 
menta p. The upper limit of the range depends on 
the number of functions N and the oscillator length b 
and, in view of the properties of oscillator functions, 
can be expressed as 


So far, to be on a safe side, we used the large ba¬ 
sis of functions {N = 300) to describe the bound 
and continuous-spectrum states. Now, we are going 
to determine a minimal set of oscillator functions, 
which gives solutions with necessary precision. We 
again turn our attention to the square-well poten¬ 
tial, because the wave function and the T-matrix for 
this potential are obtained in a simple analytic form. 
This helps us to verify the precision of our calcula¬ 


tions. The exact (391 and calculated T-matrices de¬ 


termined with N = 25, 50, and 75 oscillator func¬ 
tions are shown in Fig.[^ As is seen, 25 functions 
cannot provide with a good precision for the T-mat¬ 
rix. The calculated T-matrix noticeably deviates from 
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0 5 !0 15 20 25 30 


p, fm 

Fig. 9. T-matrix convergence for the square-well potential 
calculated for two values of energy E = 1 and E = 10 MeV 


functions is consistent with the results presented in 
Figs. [2 an di as for the convergence of calculations of 
the bound state and the phase shift. 

A similar picture is observed for the Gaussian 
and exponential potentials. Unfortunately, we do not 
know the exact T-matrix for both of them. For these 
potentials, the T-matrix calculated with N = 300 
oscillator functions can be considered as “exact”, as 
this number of functions provides us with a stable 
solution and the exact phase shift. With this defini¬ 
tion, the calculations with N = 50 and N = 75 basis 
functions are almost indistinguishable for the “exact” 
T-matrix. As for the Yukawa potential, one needs at 
least N = 125 to be close to the “exact” T-matrix. 

3.4- Asymptotics 

In this section, we consider the asymptotic behavior 
of the T-matrix as a function of n for large values 
of n. In Refs. [TJH], it was shown that the expan¬ 
sion coefficients for wave functions have the following 
asymptotic form: 

Cn = ~ (bRn). (43) 



Fig. 10. Convergence of the T-matrix for the square-well 
potential. Results are obtained for E = 20.0 MeV and with the 
oscillator length 6=1 fm 


Thus, the expansion coefficients = tnL{b, k) for the 
T-matrix can be represented as 


tn — t^iL ib-i k^ — 


V 


'i'kL) 


~ (bR„) (bR„) = U (bR„) Cn, 


(44) 


where 


Rn = V4,n + 2L + 3 


(45) 


is the turning point for a classical harmonic oscillator 
in the three-dimensional space. 

In Ref. [S] (see also EH), it was discovered that 
there is another contribution (which was called as the 
short-range (SR) contribution, while the asymptotic 
term in Eqs. (431 and (441 is called as the long-range 
(LR) contribution) to the asymptotic form, which re¬ 
lates the expansion coefficients Cn and with the 
wave function and the T-matrix 


the exact one, whereas the T-matrix calculated with 
IV = 50 and Y = 75 is sufficiently close to the exact 
T-matrix. We note that the larger the number of os¬ 
cillator functions involved in calculations, the larger 
is the range of momenta p, where the T-matrix is de¬ 
scribed by these functions. Note that this number of 


{n\^kL) ^ (-1)” V^u-^kL {Rn/b), 


(46) 


tnL= (n 


V 


^ (-1)” (Rn/b, k). (47) 


Equations (431 and ( |44| establish some relations be¬ 
tween the expansion coefficients Cn and tn and the 
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wave function and the T-matrix for the large values 
of momentum p = Rn/b. 

Let us consider the calculation of for the square- 
well potential in more details. To determine one 
has to calculate the integral 

_ a 

tn = tn (b,k) = VoAkLy '^kJdrr'^^nL{r,b) jLikor) = 

— a/b 0 

= VoAkLy Jdpp^^nL {p, 1) JL [hk^p), 

0 


where p = r/b. If a/b 1 or a/b > \/An + 2l + ?>, 
then the integral can be extended to infinity, and 
it gives the expansion coefficients for Bessel func¬ 
tions. Thus, we obtain the long-range approxima¬ 
tion. For small values of the ratio a/b 1 or a/b 
<?; \/An -\- 21 -\- we can use the approximate formula 

for oscillator functions (Laguerre polynomials). 

It was also demonstrated in jS] that, in some cases 
(which depend on the value of oscillator length b and 
the shape of a potential), the SR contribution is much 
larger than the LR one. In addition, one has to con¬ 
sider both of them in some cases. Unfortunately, we 
do not know the exact asymptotic form of the wave 
function Aikh/p) in th® momentum space and the T- 
matrix for large values of p as well. It makes difficult 
to realize the short-range approximation. 

Note that the knowledge of the asymptotics for the 
wave function and the T-matrix allowed us to for¬ 
mulate strategies (see Ref. jj and also Ref. uni) to 
obtain convergent results for the phase shift with a 
minimal set of oscillator functions. 

We have explicit forms of the wave function and 
the T-matrix only for the square-well potential. We 
will use it to check the asymptotic behavior of the 
T-matrix. 

In Fig. 11 we compare the calculated (E) and 
asymptotic long-range (A(LR)) forms of the T-matrix 
with the exponential potential. As is seen, the asymp¬ 
totic long-range form is valid for small values of 
b. Starting from b = 1.0, the asymptotic form is much 
smaller than the exact form. A similar picture is ob¬ 
served for the Gaussian potential. However, for the 
Yukawa and square-well potentials, the long-range 
form gives a very small contribution comparing with 
the exact form. 

In Fig. we demonstrate how the asymptotic 
short-range approximation works for the square-well 



n 

Fig. 11. Exact (E) and asymptotic (A(LR)) forms of expan¬ 
sion coefficients tn calculated with the exponential potential 
for the energy E = 10.0 


potential. We show results for two values of oscillator 
length b: the smallest (6 = 0.5) and largest (6 = 2.0) 
ones. As is seen, the asymptotic short-range form is 
valid for small values of b. For large values of b, the 
asymptotic form coincides with the exact one for the 
whole range of the quantum number n. 

It should be stressed that the long-range asymp¬ 
totic form gives zero contribution to the presented 
results due to a specific shape of the square-well po¬ 
tential. 

4. Conclusion 

We have studied the properties of the T-matrix in the 
discrete oscillator representation. It is demonstrated 
that the T-matrix in the oscillator representation can 
be presented in the vector and matrix forms. The vec¬ 
tor form is suitable for investigating the T-matrix on 
the half-on-shell space, while the matrix form is more 
appropriate for the full off-shell space. The set of lin¬ 
ear equations for the T-matrix expansion coefficients 
is deduced for both vector and matrix forms. 

We calculated the T-matrix expansion coefficients 
for four different potentials. It is shown that the 
rate of decreasing of the T-matrix expansion coeffi¬ 
cients depends on the shape of a potential and on 
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Fig. 12. Expansion coefficient for the exact T-matrix (E) 
compared with the asymptotic short-range form (A(SR)). Cal¬ 
culations are made for E = 10.0 with the square-well potential 

the oscillator length. It is also shown that the T- 
matrix expansion coefficients slightly depend on the 
continuous-spectrum state energy. We recall that the 
energy of scattering states is considered in the range 
0 < i? < 20, which is typical of the low-energy nuclear 
processes. 

It is shown that the calculations of the T-matrix 
in the discrete representation is a reliable way for ob¬ 
taining the information on the behavior of a quantum 
mechanical system. 

We have investigated thoroughly the asymptotic 
properties of the T-matrix in the discrete space and 
established a relation with its asymptotic form in the 
continuous coordinate and momentum spaces. 

The present discretization method can be easily ex¬ 
tended for the T-matrix of real physical systems such 
as many-channel and many-cluster systems. 
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B. C. BacuAeecuKuu, M.ff. CoAoxa-KAUMuaK 

T-MATPHIJH B AMGKPETHOMY, 

OGIJHJIHTOPHOMy HPEACTABJIEHHI 

P e 3 K) M e 

Jl^ocjiifl:»ceHO BjiacTHBOCTi T-MaxpHpi /i;jih SB’nsanHx CTaniB Ta 
CTaniB nenepepBHoro cneKxpa y ^ncKpeTHOMy, ocu,HjmTopHO- 
My npeflCxaBjieHHi. ^ocjii;],}«;eHHH npoBO^^BTbCH flJiH MO^ejib- 
Ho'i npoSjieMH — uacTHHKa b nojii penTpajibHoro noTeHii;iajiy. 
BHBe/i;eHO cncxeMy jiininHHx piBHHHb, P03b’h30k nxnx BH3Ha- 
uae Koe4>iLi,ieHTH po3Kjia;],y T-MaTpHii;i no ocpHJiHTopHHx 4)yH- 
Kpinx. Mh BubpajiH uoTupn noTenpiajin (rayciBCbKnn, loxa- 
BiBCbKHH, eKcnoneHpiajibHHH xa noxenpiaji npnMOKyxnoi hmh) 
fljin /i;eMOHCxpai];ii ocobjiHBOCxen T-Maxpnpi xa i’l aajiejKuocxi 
Bi/i; 4>opMH noxeHii.iajiy. Mn xaKO»c BHBuaeMO hk Koe4>iu;ieHXH 
po3Kjia/iy T-Maxpnpi 3ajie:»caxb Bi;], napaivrexpiB ocpnjiHxopno- 
ro 6a3Hcy — ocpnjiHxopHoi' ^ob»:hhh xa uncjia 6a3HCHHx 4)yH- 
Kpin, aajiyneHHx y po3paxyHKax. 
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